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ABSTRACT 



Although still poorly understood, the chemistry that occurs on the surfaces 
of interstellar dust particles profoundly affects the growth of molecules in the 
interstellar medium. An important set of surface reactions produces icy mantles 
of many monolayers in cold and dense regions. The monolayers are dominated 
by water ice, but also contain CO, CO2, and occasionally methanol as well as 
minor constituents. In this paper, the rate of production of water-ice dominated 
mantles is calculated for different physical conditions of interstellar clouds and for 
the first time images of the morphology of interstellar ices are presented. For this 
purpose, the continuous-time random-walk Monte Carlo simulation technique has 
been used. The visual extinction, density, and gas and grain temperatures are 
varied. It is shown that our stochastic approach can reproduce the important 
observation that ice mantles only grow in the denser regions. 

Subject headings: ISM: abundances-ISM: molecules-dust-stars:formation 



1. INTRODUCTION 



In cold and dense re gions of interste 



istence of m antles of ices (|Williams et al 



constituent ( Whittet et al 



1998 



lar cl ouds, infrared absorption studies show the ex- 
1992), of which water ice is observed to be the main 



Pontoppidan et al.l 12004 ). The formation and destruction 



of the interstellar ices is particularly important in regions where new stars and planets are 
born. During the first stage of star formation, small cores in a cloud assembly collapse under 
their own gravitational force, resulting in high gas densities and low temperatures. Under 
such conditions, virtually all gaseous species heavier than hydrogen and helium accrete onto 
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the grains, where they can undergo chemical processes. I n protostellar regions, the ices evap- 



orate and change the nature of the gas-phase chemistry (ICaselli et al.lll993t iBottineUi et al. 
2004h . 



The formation of water ice-dominated mantles is not yet fully understood. Recombina- 
tion of smaller species on the surface and condensation of complete molecules from the gas 
phase are the competing processes, but quantitative estimates of their efficiency under inter- 
stellar conditions are lacking. Also, the reason for the apparent absence of water ice in more 
diffuse clouds, which are exposed to ultraviolet and visible photons, is still not clear. From 
observations of these regions, we know that the ice mantles constitute approximately less 
than a few monolayers of water, the cur rent detection li mit. A recent attempt to understand 
the formation of ice layers was made by lPapoularl (120051 ) who concluded that direct accretion 
of water from the gas is necessary to initiate ice-mantle formation. His main explanations 
for the lack of water ice in the diffuse regions are the high grain temperature and low flux 
in these regions. 



As iPapoularl (120051 ). our goal in this paper is to explain the growth of monolayers of 
ice in dense sources and the absence of this growth in diffuse sources. With the use of a 
so-called microscopic Monte Carlo approach (see below), we focus on the accretive recombi- 
nation mechanism using a set of surface reactions undergone by species that either accrete 
onto a grain surface or are products of other surface reactions that remain on the grain. 
The reactants are limited to species containing the elements O and H. We also consider 
photodissociation processes for surface species caused by ultra-violet photons. In diffuse 
sources, these photons are mainly those of the external radiation field, whereas in dense 
sources, the much smaller flux of photons arises indirectly from cosmic ray bombardment of 
H 2 , w hich produces ions and electrons. The electrons can then excite H 2 , leading to fluores- 
cence ( jPrasad fc Tarafdarlll983l ; iRuffle fc Herbstll200ll ). In translucent and dense clouds, we 
find that ice mantles grow and we follow their morphology in some detail. To the best of our 
knowledge, this is an initial attempt to explain the morphology of ice grown via chemical 
reactions. 



2. MONTE CARLO METHOD 



The exact details of the continuous-t ime random-walk ( CTRW) Monte Carlo m ethod 
are already explained in previous papers (jChang et al.ll2005l ; ICuppen fc Herbstl 120051 ) . We 
therefore give only a brief summary here. The initial grain surface is divided up into a 
square lattice of adsorption sites, each of which corresponds to a local energy minimum. We 
start from a bare surface of amorphous carbon with a high degree of roughness. Carbon 
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is used rather than olivine because of the more severe lack of available desorption data for 
the adsorbates we must consider on the latter surface material. Since the major portion of 
carbonaceous grains is of an amorphous nature rather than crystalline, we treat the grain as 
being amorphous. This choice also means that the grains do not have a metallic character, 
as they would had they consisted of graphite. Figure [1] shows the topology of the surface of 
the bare carbonac e ous g rain we use throug hout the paper . It is t he same as surface (d) in 
Cuppen fc Herbstl (120051 ) and surface III in ICuppen et al.l (120061 ). ICuppen fc Herbstl (120051 ) 
explain the procedure used to obtain this surface. The roughness is caused by a variety of 
surface imperfections including vacan cies (holes in the su r face), edges, islands, and ad-atoms, 
all of which are shown in Figure 2 of ICuppen &: Herbstl (120051 ). The species are confined to 
a square lattice structure even though we intend to simulate the growth of an amorphous 
ice mantle. An off-lattice simulation method, where the species are not confined to regular 
lattice positions, would seem to be a better choice, but this procedure is computationally 
expensive. Since we do not have a potential or force field to reliably describe the system, 
it would not result in much more realistic simulations. The surface plotted in Figured] has 
50 x 50 lattice site s. Because this size is above the limit where finite size effects play a role 
(jChang et al.l 120051 ) . scaling the results to larger sizes does not cause a problem. Indeed, we 
scale the results to a standard grain with a radius of 0.1 /im, which corresponds to a 250-500 
times larger grid, depending on the site density chosen. The smaller grid was used to keep 
the cpu time for the simulations reasonable. 

Atoms of a species labeled A are deposit ed on these grains during our simulations with 
a deposition rate in monolayers (ML) s _1 of (IBiham et al.ll200ll ) 
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where n\ is the absolute gas abundance of species A, p is the surface site density, and v is 
the mean velocity of the atoms and molecules in the gas: 
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8kT, 



gas 



nrriA 



(2) 



with T gas the temperature of the gas and rriA, the atomic mass. We use p = 5 x 10 14 cm -2 
in diffuse clouds and p — 1 x 10 15 cm" 2 in dense clouds. The latter value is obtained from 
measurements of high density ice ( Jenniskens et all 1995 ). The value used for diffuse clouds 



was determined for amorphous carbon (IBiham et al. 



200ll ). In our simulation, atoms hit the 



surface under a randomly determined angle. The atom will stick to some environment on 
the rough surface depending on its incoming angle. Once on the surface, the atom can hop 
over the surface, or desorb. The rates (s _1 ) of hopping and desorption are respectively 
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and 



R P 



v D exp - 



E 



D 



(4) 



where T is the surface temperature of the grain, Eb and Ejj are the hopping barrier and 
desorption energy (K), and v h and v D are the attempt frequencies for hopping and desorption, 



respectively. The frequencies Vb and vr> are both taken to be 10 
value for physisorbed species (IBiham et al.ll200ll ). 



12 



which is a standard 



2.1. Binding, Hopping, and Activation Energies 



We assume that species close to step edges or other protrusions of the surface are more 
strongly bound because of the extra lateral physisorption "bonds" with the surface. Only 
lateral interactions with the water ice and the bare carbonaceous grain are considered. 

For the binding energy of a particular adsorbate species, labeled A, we use the equation 



i K2 o) = £ S A + ai c E A + az H2 o^H 2 o> 



(5) 



where E^ is the binding energy of A with the particular surface lying beneath it, E^ is the 
binding energy of A with amorphous carbon, E# Q * s the binding energy of A with water-ice, 
a is the fraction of the binding energy that characterizes lateral binding, and i s (s = c, H2O) 
is the number of horizontal neighbors of the particle with carbon or ice. For the hopping 
energies, we use the equation: 



E b («c, i H2 o) 



0.5E A + ai c E A + ai n2 oE^ 20 , 



(6) 



which is similar to equation (jSj) except for the factor 0.5 in the first term. Very little is known 
about the relation between the hopping barrier and the desorption energy. For atomic hydro- 
gen on amorphous carbon and polycrystalline sil icate surfaces a fa ctor between the hopping 
barrier and desorption energy of 0.78 was found (IKatz et al.lll999l ). This value was obtained 
by fitting rate equations to temperature-programmed-desorption ( TPD) experiments. In 
most astrochemical models, a much lower value of 0.3 is assumed ( ITielens fc Allamandola 
19871 ). Experiments of water-on- water diffusion indicate a similar value to .3 whereas CO 



on water-ice experiments suggest a slightly higher ratio (jCollings et a 



2003). Our value of 



Ulbricht et all f|2002 h 



0.5 lies closer to this last study and corresponds to the value found by 
for oxygen adsorption on carbon nanotubes. The higher value found by IKatz et al.l (jl999l ) 
could be the result of an average over different hopping barriers. Because our method of 
lateral bonds already accounts for this effect, we take a lower value. Appendix IA1 discusses 
the effect of this hopping factor on the results. Note that the lateral bonds must be broken 
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in a hopping movement and that we do not account for the bonds that are re-formed since 
we assume the barrier to be independent of these bonds. We do not consider a possible 
vertical interaction upwards for embedded layers. If a particle has "upper neighbors" we do 
not allow it to evaporate but it can hop, leaving a vacancy, or "pore" , in the material. 

The a factor is taken to be 0.2 in all simulations. Once again, very little is known 
about what this factor should be for physisorption. We choose our relatively low value to 
be conservative; non-zero values are needed to explain the high efficiency of H 2 formation 
at relevant surface temperatures in diffuse clou ds but high values can l ead to complex TPD 
spectra in the laboratory with multiple peaks (ICuppen fc Herbstll2005l ). The factor has the 
largest effect on the hopping and desorption of hydrogen and on the hopping of oxygen atoms. 
All other species have such high energies for hopping and desorption that they remain in 
their sites independently of any changes in energies due to a different a factor. 

If two particles encounter each other due to hopping and land in the same lattice site, 
a reaction can occur. Table [1] lists the surface reactions used to study ice formation, their 
estimated activation energies, and a factor j3 that stands for the fraction of products remain- 
ing on the surface after a reaction. All of the surface reactions in Table [TJ are exothermic. 
The excess energy that is released in an association reaction (a process with one product) 
can cause the product to desorb from the surface. To the best of our knowledge, (3 has been 
determined e xperimentally on ly for the formation of H2 on olivine (0.33) and amorphous car- 
bon (0.413) (IKatz et al.lll999l ). but not on ice. The (3 value for the formation of water from 
H and OH can be extracted from molecular dynamics simulations of the photodissociation 
of solid water performed by iKroes fc Anderssonl (120061 ). They find that a fraction of the H 
+ OH photodissociation products reacts back to water, 0.9 % of which is released into the 
gas phase. We take this value of (3 = 0.991 for all single-product association reactions. Only 
products in the outer layers, which are directly in contact with the gas, a r e allowed to evap- 
orate. A value of (3 near unity is also found in the theory of iGarrod et al.l (120061 . 120071 ). This 
fraction is mainly important to remove water from the surface in diffuse clouds. It is of less 
importance for the H2 production, because H2 can evaporate thermally. Since the binding 
of molecular hydrogen to amorphous ices is similar to its binding to amorphous carbon, the 
/?-value is expected to be similar on this substrate. We will use the value of (3 — 0.413 for 
this reaction. A high er value will result in a slightly lower hydrogen recombination efficiency 
JCuppen et alil2006h . 



In contrast with might be expected from the case of gas-phase reactions, surface reactions 
with activation barriers can have noticeable contributions in our treatment, because the 
number of reaction trials the species can get before one of them hops away multiplied by the 
probability of reaction can reach unity (lAwad et al.ll2005l ). In our approach, if two species 
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land in the same box and can react with activation energy, we do not follow the details of 
the competition but determine immediately which process - reaction or hopping - occurs 
by comparison of the two rates. There is experi mental and/or theoretical evidence that 



these reactions do indeed occur on cold surfaces (IKroes &: Anderssonl 120061 ; iHiraoka et al. 



1994 ; Tielens fc Hagen 1982 ; Watanabe fc Kouchi 20021 ). Here we assume that the reaction 



barrier can be treated as a rectangular potential barrier, through which tunneling occurs 
with a rate (s -1 ) 



R 



react 



U r P, 



react 



exp(-^ v / (2^ a )) 



(7) 



where a is the width of the barrier, /i the reduced mass, v r the attempt frequency for reac- 
tion, and E a the reaction activation energy. We assume that the trial frequency for hopping, 
desorption and reaction is the same. We take here a width of 1 A. The width of the dif- 
fusion barrier is probably larger (D. Woon, private communication) and diffusio n therefore 



most likely occurs through thermal hopping, as in the case of H 2 formation (IKatz et al. 



19991 ). The difference in treatment between diffusion and hopping has a large influence on 



the temperature dependence of reactions with a barrier. The probability P re act is temper- 
ature independent according to equation (J7|) whereas the thermal hopping rate, which is 
inversely proportional to the residence time, increases with increasing temperature. Tunnel- 
ing leading to reaction therefore becomes more important as the temperature is lowered, but 
the decreased hopping rate makes it harder for two reactants to find each other. Detailed 
experiments and calculations need to be carried out to give more information about how to 
treat this type of surface reaction. 



2.2. Long-Range Diffusion and Photodissociation 

If a species tries to hop to a neighboring site that is already occupied and the two species 
do not react, we allow the particle X to diffuse to the next available free lattice site in the 
same layer with a probability P given by the empirical equation 

p = ex P (- £c(x 2 ;; H20) a, 

where Ax is the distance in sites to the next available free site, and Ed is the binding 
energy in K. This equation is chosen so that small, weakly-bound species such as atomic and 
molecular hydrogen can diffuse thro ugh porous ice str u cture s, with a probability that rapidly 



decreases with increasing distance. lAndersson et al.l (120061 ) showed in molecular dynamics 



simulations that atomic hydrogen can move through the ice matrix after photodissociation of 
water, when the hydrogen atom is translationally excited. Eq. [8] allows this motion to occur. 
The probability is an estimate based on the arguments mentioned above, since information 
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for thermalized hydrogen at in the range of 10 to 20 Kelvin is lacking. Eq. [8] does not 
contain an attempt frequency because this probability is only called if a normal hopping 
event cannot occur since the neighboring site is already occupied. The attempt frequency 
is therefore already in the hopping rate as is the temperature dependence. The particle can 
diffuse in all six directions. If the particle does not diffuse, it remains in its original position. 
Using binding energies discussed in the next section, the probability of diffusion of H over 
two sites P is 4.8 x 10~ 3 while it is 1.5 x 10~ 8 for the two-site diffusion of a heavier ozone 
molecule. 

Surface species can be destroyed by photodissociation, either directly or by cosmic-ray- 
induced photons. We assume that the photodissociation rate of the surface species is the 
same as in the gas phase, and is given by the standard expressions 

R = aphotoexpi-JphotoAv) (9) 

and 

R = a cr ( (10) 

for direct and cosmic-ray-induced photodissociation, respectively. The constants are given 
in Tableland for the cosmic ray ionization rate £ we use 1.3 x 10~ 17 s -1 . These assump- 
tions are only parti a lly in agreement with molecular dynamics calculations for ice done by 



Kroes fc Anderssonl (120061 ). who show that the photodissociation cross section for water ice 
does not have exactly the same dependence on photon energy as is the case for gas-phase 
water. In our simulation, all molec ules in the grain man tle can be photodissociated and a 



shielding mechanism is not included. iNguyen et al.l (120021 ) assumed a limit of 100 monolayers 



that could be reached by photodissociation. In our simulations, only in dense regions where 
photodissociation plays a minor role can these thicknesses be reached. 

We further assume that the heavy product fragments remain on the surface after pho- 
todissociation. These product fragments ca n, however, react leading t o the possibility of 



desorption in association reactions. Following [Kroes fc Anderssonl (120061 ) we assume that 40 
% of the atomic hydrogen formed by photodissociation of H2O evaporates. This value is an 
average of the value they obtained for the first three layers. We also include their indirect 
mechanism for the photodesorption of water in which a thermally excited hydrogen atom 
kicks out one of the water molecules with an overall probability of 0.04%. 



3. ENERGIES 



Table [3] contains the binding (desorption) energies used between assorted adsorbates and 
surfaces in our calculations. Many interaction energies are needed to describe the system, 
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not just between the different species with the bare carbonaceous substrate but also with 
ice and with the other adsorbed species. The problem is, however, that many of these 
interaction energies are poorly known. For amorphous carbon, as far as we know, only 



the desorption en ergies for H and H 2 are experimentally determined (jPirronello et al.lll999 
Katz et al.l Il999l ). with values of 658 K for H and 542 K for H2. Additional studies 



have been undertaken with graphite and carbon nanotubes. Even though our model grain 
consists of a morphous carbon, we occasionally have to rely on studies based on other forms 



of carbon. iGhio et al.l (Il980l ) calculated the van der Waals interaction between graphite 



and atomic hydrogen and found an energy of 503 ± 6 K, which i s in excellent agreement with 
the overview of experimental data given by IVidali et al.l (ll99ll ). in which the best estimate 
is given as 499 ± 3 K. For th e H2 i nteractions with carbonaceous material we found two 
additional studies. lHan &: Led (120041 ) used DFT method s to ca lculate a value of 400 K for 
H 2 on carbon nanotubes, while the survey of IVidali et al.l (Il99ll ) gives a best estimate of 660 
± 6 K for H2-graphite. In both cases, we use the amorphous carbon results. 

For the other species, there are no amorphous carbon data and we have to make use of 
theoretical or ex perimental results for graphite, carbon nanotubes, or other forms of carbon. 
Papoularl (120051 ) performed chemical simulations to obtain the interactions between O, OH 
and H 2 and typical hydrocarbon surface functional groups. He found energies of 300 K and 
500 K for O and OH, respectively, while for water he found a range of energies between 250 
and 1000 K. The O and OH values seem rather low compared with the energies for the much 
lighter species H and H 2 . The calculations were performed between small single molecules 
and the absorbate. In reality, there will probably be many more contributions from other 
atoms in the substrate resulti ng in a higher inte raction. We prefer to us e the estimates of 
800 K jTielens fc Hagenlll982h for O and 1360 K (Urlen fc Robinsonl 1 1975^ for OH. For H 2 0, 
there is another theoretical value by Picaud et al. ( 2004T). who d etermined th e energy to be 
800 K . Experimental results give 1870 K (Uvgul & Kiselev!ll97oh and 2150 K (bale & Beebe 
1964 ) for the same s ystem. We use 2000 K based on the experiments. For 2 , we use the 
value of 1440 K from lUlbricht et al.l (120021 ) obtained by TPD experiments on graphite, while 
for the species O3, 2 H, and H 2 2 , we simply add the binding energies of their constituent 
parts: 



E D {o 3 ) = E D {o 2 ) +E D (p), 
E D (o 2 n) = E D (on) + E d (o,) 
E D (n 2 o 2 ) = E d (o 2 h) + E d (h). 



(11) 
(12) 

(13) 



Let us next consider the energies of the species on solid water/ice. The problem here is 
that there are many types of ice, both in structure and in morphology and porosity, which 
makes it difficult to compare different experimental results. The porosity has two effects: 
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(i) the particles in the pores cannot desorb directly but will perform a random walk on 
the wall of the pores until they reach the outer surface, and (ii) the increased roughness 
in the pores means that the number of binding sites with high energy is probably larger. 
Since we already account for the latter effect by considering multiple types of sites, we 
use non-porous resul ts as much as possible . A considerable amoun t of data is available for 
molecular hydrogen. iHornekaer et al.l (120051 ) and lDulieu et al.l (120051 ) find a large distribution 
of energies using TPD experi ments. If we only consid er the mean value for the non-porous 
substrate we get 4 40 K from IHornekaer et al.l (120051 ) for D 2 and 520 K and 550 K from 



Dulieu et al.l (120051 



value obtained by 



for H 2 and D2, respe c tively. These values are close to the theoretical 



Holl enbach fc Safpeterl (119701 ) of 550 K and the experimental value of 



Sandford fc Allamandolal (19931) of 550 ± 35 K, based on a mixture of water and methanol 
ice. The results of IPerets et al.l ( 20051 ). however, are in disagreement with this cluster of 
values. Performing TPD experiments for HD and D 2 on high density ice (HDI) and low 
density ice (LDI), they found three types of sites for the LDI, which has a higher porosity 
and one type of site for the HDI with binding energies 798 K (HD) and 836 K (D 2 ). We 



use 44 K, the lower limit of the distribution and the peak energy found by IHornekaer et al. 



( 120051 ). Since we include higher energy sites near step edges in our model, we still obtain a 
distribution of desorption energies in this way. 



F or H on ice, we found t h ree th eoretical values: iHollenbach &: Salpeterl (jl970l ). lBuch &: Czerminski 
dl99lh and kl-Halabi etaD (boQ2h determined energies of 45 K, » 500 K and 400 ± 50 



K, respectively. The values of 



Ho 



are for crystalline ice whereas the 



lenbach fc Salpeter 



( 1970 ) and lAl-Halabi et all (12002( 1 



Buch fc Czerminskil (Il99ll ) value is for amorphous ice. 



Perets et al.l (120051 ) determined a significantly higher val ue, 720 K, from TPD exper i ments 
on HDI. We use 650 K, which is a new value calculated by lAl-Halabi fc van Dishoecki (120061 ) 
for the binding energy of atomic hydrogen on amorphous ice. They found that the binding 
on amorphous ice is significantly higher than on crystalline ice due to the surface structure. 
Even a surface without protrusions of amorphous ice has several small cavities that will bind 
the hydrogen more strongly whereas the crystalline ice is even flat on an atomic scale. 



For water on amorphous ice we use 5640 K, which is the result obtained bv lSpeedv et al. 



( 1996 ). and is in excellent agreement with a value of « 5600 K found by iFraser et al. 
(l200lf). For Q 2 we us e 100 K, which is an ave rage value obtained from the TPD data 



by lAvotte et al.l (120011) and Collings et al.l (120041) . For O and O H we again take the esti- 
mates from lTielens fc Hagenl (119821 ) and lAllen fc Robinson! (119751 ). respectively, and for O3, 



2 H, and H 2 2 , we again add the binding energies of their constituent parts. 

Let us now cons ider H 2 as a substra te, since, as we will see, its abundance in the mantles 
can be considerable. IVidali et al.l (Il99ll ) give an estimate of the desorption energy for H from 
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H2 of 35.4 ± 5.8 K + kT which leads to 45.4 K at 10 K. The rest of the desorption energies 
from H 2 are, starting from the ice desorption energies, scaled according to the equation: 

„ , v En(X on HoO) „ . . 

E D (XonR 2 ) = ^£ D HonH 2 . 14 

£/£)(H on H 2 Oj 

Once again, energies involving O3, O2H and H2O2 are formed from their fragments. 

The reason the desorption energy of O2H is considered to be the sum of the desorption 
energies of OH and O instead of O2 and H, which are better known, is the capability of O2H 
to form H-type bonds. We are aware that desorption energies obtained in this manner have 
a large uncertainty. Fortunately, our simulations show the uncertainties in the desorption 
energies of the strongly bound species to have very little effect on the final result, since 
these species will remain on the surface at these temperatures within the time-scale of the 
simulation regardless of their exact desorption energies. 



4. RESULTS 

In all simulations the water was formed through surface reactions and not by direct 
accretion from the gas phase, even though gas-phase water can be detected in some dense 
objects. This exclusion is easily justified. A typical value obtained in gas-phase models for 
the fractional abundance of water is at most 10~ 6 . For a cloud of density 10 4 cm -3 , and a 
standard site density, we then calculate that the flux rate for gaseous water at 10 K upon 
a grain is 3 x 10~ 14 ML s -1 , so that it takes roughly 10 6 yr to deposit a single monolayer 
and the enormous time of 10 8 yr to deposit the 100 monolayers typically associated with 
dense cloud mantles. The time scales are much longer for diffuse clouds. The flux of oxygen 
atoms arriving on grains is much higher than the flux of water molecules, typically by a 
factor of 100 in dense clouds. With a high efficiency of conversion to water by surface 
reactions, 100 monolayers of water can then be produced in 10 6 yr, in reasonable agreement 
with our dense-cloud results discussed below. So, achieving significant mantle abundances 
in reasonable times does require surface reactions. 



4.1. Diffuse and Translucent Clouds 

We have performed a set of simulations for different values of the visual extinction Ay, 
total hydrogen density (wh), and gas and grain temperatures. Figure [2] plots the surface 
abundances in monolayers for the different species in our model as functions of time. In this 
figure, the results for diffuse and translucent clouds are contained. All panels are labeled with 
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a letter followed by a number. The panels starting with the same letter are grouped in rows 
and have the same visual extinction, which increases from top to bottom. The digits along 
a row indicate the different physical conditions. For the three leftward panels, the density is 
fixed but the temperatures decrease as one goes from left to right, while the rightmost panel 
(e.g. A4) has the same temperature as the second of the three left panels (e.g. A2) but a 
higher density. The exact physical conditions for each of the simulations are given in Table 
HJ In all cases shown in Figure [2j we assume the clo ud gas to be initial ly atomic, with an 



oxygen abundance of 3.0 x 10~ 4 n# that of hydrogen (jNguyen et al.l 120021 ). The calculations 
are run for sufficiently short times that the gas-phase abundances remain unchanged. In all 
cases we start with a bare carbonaceous grain. 

Figure [2] clearly shows that in diffuse regions (Ay = 0.5, 1, rows A and B) less than 1 
ML of water-ice is present at a time of 3 x 10 5 yr, the largest time studied. For panels A1-A3, 
the water-ice abundance has reached steady-state conditions while in panel A4, which refers 
to a higher density, steady state seems to be reached around the end of the simulation. In 
panels Bl and B2 the water layer is still growing, albeit slowly, while in panels B3 and B4 
the water abundance is decreasing. Although water-ice is not destroyed chemically in our 
model, a steady-state results when destruction by photons and desorption from the mantle 
balance formation. Water-ice is the most abundant species for the diffuse clouds in rows A 
and B; nonetheless, other species can be seen in the panels. Despite its low desorption energy, 
molecular hydrogen can reach mantle abundances of a significant fraction of a monolayer in 
diffuse clouds. In the translucent clouds, rows C and D, the water-ice abundance exceeds 
1 ML by the end of 3 x 10 5 yr and is still growing. Water-ice is dominant in all panels 
except D4, where OH has a comparable total abundance. As far as we know, no data from 
observations is available about the abundance of OH present on grains under comparable 
conditions. 

If we look at the influence of the physical conditions on the formation of ice for diffuse 
and translucent clouds, it can be concluded that the gas and grain temperature affect the 
ice formation mechanism only to a small extent. In particular, a low surface temperature 
results in a slightly thicker ice mantle. In general, the ice thickness is mostly determined by 
photodissociation under the influence of the radiation field and by the gas density. Thus, 
as one goes down a column and to the rightmost panel of any row (higher density), the 
water-ice abundance increases. Figure [2] shows a clear difference in the composition of the 
ice layer depending on the visual extinction. At Ay = 0.5 (A panels) the mantle mainly 
consists of water, and molecular hydrogen. For increasing Ay the amounts of OH and O3 
increase. The surface abundance of atomic oxygen appears to peak at Ay = 2 (C panels). 

Fig. [3] shows a vertical slice of the ice at 3 x 10 5 yr for four different conditions - panels 
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A2, B2, C2, and D2 - in which visual extinction and gas density increase. The different 
molecules are color coded in the panels, while gray represents the carbon substrate and 
black represents the void. One sees immediately that the mantle for A2 is virtually non- 
existent, while that for D2 has grown to a significant number of monolayers. The mantles 
are certainly porous; there is a lot of black interspersed among the color-coded molecules. 
Due to the strong radiation in all simulations, the ice goes through many cycles of formation 
and destruction. Let us take a more detailed look at mantles C2 and D2. The top layers 
of the mantles are primarily water (white color). As one goes deeper into the mantle, 
one notices more green (OH), and finally towards the bottom one finds lots of magenta 
(ozone). In understanding this dependence of abundance on mantle height, the topology of 
the mantle is important. Blocking of the evaporation of the hydroxyl radicals and oxygen 
atoms after photodissociation and obstruction of the deposition of fresh hydrogen from the 
gas phase are the main reasons that allow the large abundances of OH and O3 to be formed 
away from the surface. Atomic hydrogen leaves the surface either by direct desorption after 
photodissociation or in the form of H2 that can diffuse through the pores. We used the same 
desorption probability of H after water photodissociation of 0.4 throughout the ice. If H 
atoms accreting from the gas have a greater penetration depth or if a depth dependence of the 
desorption probability upon photodissociation was considered, the large abundance variance 
with depth we predict could be partially reduced. Unfortunately, very little experimental 
information on the structure of water ices formed through surface reactions under interstellar 
conditions is available so that our images of the ice mantles cannot be easily compared with 
real ices. 



The common gas-grain codes (IRuffle fc Herbstll200ll ) do not keep track of the molecular 



structure of the ice layers so structural variations cannot be found. The mean field approach 
of most of these methods allows all species to react with each other independently of their 
position. The mantle abundances of the hydroxyl radical and ozone will therefore be lower 
in these simulations using the same surface reaction network since their reactions with H 
atoms will occur more readily. Even in our variegated ices, the mantle abundances of OH 
and O3 may decrease to some extent if other reactive species like CO are included and are 
allowed to react with OH and O3 in our reaction network. For the latter case these reactive 
species must be able to diffuse in the lower layers of the ice. 

Molecular hydrogen formation from the recombination of two hydrogen atoms on a 
surface is the most important reaction in the diffuse interstellar medium. The recombination 
efficiency for H 2 formation is defined by the ratio of twice the hydrogen molecules formed 
and leaving the surface divided by the flux of hydrogen atoms. This efficiency has been 
determined over the time interval that steady state is achieved for the 16 sets of diffuse and 
translucent physical conditions (A-D) given in Table H] and the results are given in Table 
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[5] by the left value in each column. The values on the right indicate the efficiency for the 
same densities and temperatures but without oxygen in the system, which means without 
ice being formed. There is a clear difference in efficiency, with the ice surface being less 
efficient in most cases. This is due to two effects: the first is that the hydrogen atoms are 
also used to form water, but this is only of minor importance. The larger effect is caused 
by the changing energetics between the water layers and the carbonaceous substrate. In 
particular, the binding energy of H on an ice surface is slightly less than on a carbonaceous 
grain and also the surface roughness changes due to the build-up of the water layers. In these 
calculations, we only use an a of 0.2 for the lateral bond strength. A higher a would result 
in a more efficient molecular hydrogen production because it ca uses the hydrogen atom s 
to be more strongly bound and to reside longer on the surface (ICuppen fc Herbstl 120051 ) . 
The efficiencies for lower surface temp eratures (12-15 K) are generally still high enough 
to produce a reasonable amount of H 2 ( ICuppen et al.ll2006l ) but the efficiencies are clearly 
degraded at higher temperatures. For smaller grains, we expect stochastic heating to reduce 
the formation of ice since both oxygen and hydrogen atoms will usually evaporate when 
a photon hits a grain. This effect wou ld cause the recombi nation efficiency for molecular 
hydrogen to approach those reported in ICuppen et al.l (120061 ) . 



4.2. Dense Molecular Clouds 

Figure H] is similar to Figure [2] but shows results for higher visual extinctions (Ay = 
5, 10) and densities (5.0 x 10 3 cm -3 - 5.0 x 10 4 cm -3 ). The maximum time considered is 
1 x 10 5 yr, at which time steady-state abundances have not yet been reached for the major 
constituents. Again the visual extinction is varied vertically while temperature and density 
for a fixed extinction are varied horizontally. The exact conditions are once again given in 
Table HI We assume that the hydrogen in the gas phase is mostly molecular. A relative 
abundance of 1 x 10~ 4 n# for atomic H is used. Although it is assumed that the mantles 
grow upon a carbonaceous substrate, it is more likely that they grow upon ice mantles of 
the type discussed in the previous section for translucent clouds. It can be seen in Figure @] 
that in seven out of the eight panels, the growth of water-ice dominates the mantles, while 
in one panel (F4) the abundance of H2 overtakes that of water at later times. Mantles of up 
to 75-80 monolayers are developed. 

In the panels, a general tendency can be observed that a lower temperature (compare 
El, E2, and E3) and a lower density (compare E2 and E4) result in less ice. This is in 
contrast with the results in Figure [21 where the amount of ice increases for decreasing surface 
temperature. In the diffuse clouds the grain temperature is in the regime where H atoms 
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easily evaporate. Reducing the surface temperature will therefore result in more hydrogen 
on the surface increasing the efficiency of the hydrogenation reactions. The lower grain 
temperature in the dense clouds lies in the regime where the desorption of hydrogen molecules 
starts to slow down and these molecules occupy a large part of the mantle, slowing down 
the rate of hydrogenation reactions. Although laboratory simulations of interstellar ice 
reactions show that hydrogenation reactio ns are still possib le in the presence of H2 because 



the hydrogen can still penetrate H 2 layers (IFuchs et al.ll2007l ). the reactions are slowed down 



which results in an optim um temperature fo r the reaction rate. For CO + H, this optimum 



lies between 12 and 15 K (IFuchs et al.l 120071 ) 



Since the atomic hydrogen abundance is much lower than the molecular hydrogen abun- 
dance, the main route to form water ice changes in dense clouds. Table E] gives the contri- 
butions of the three reactions that form water: H + OH, H + H2O2 and OH + H2. This 
table clearly indicates that the main route for formation is through the association of H and 
OH in the diffuse and translucent clouds. Since H is abundant in the gas here and since this 
reaction is barrierless, this mechanism is the most apparent one. However in dense clouds, 
the reaction between H2 and OH becomes the most common route and also the reaction of 
H with H2O2 makes a significant contribution. It seems surprising that in a dense molec- 
ular cloud a reaction with an activation barrier is so important for the formation of water 
ice. This is different fro m the gas-grain netwo rk simulations where the barrierless reaction 



remains the main route ( iRuffle fc Herbstll200ll ) along with some other surface reactions that 



were not included in this work. The difference is mainly due to the difference in treatment 
of the activation barriers, as explained earlier: because of our approach, in which the com- 
petition with diffusion is modeled directly, reactions with a barrier can have enhanced rates. 
In the approach used in the models, tunneling under the activation energy barrier is simply 
treated as a factor in the expression for the rate coefficient. 

Figure [5] gives two vertical snapshots of ice mantles. Again the surface structure is 
very rough and porous. Unlike the diffuse cloud ices, however, the structures demonstrate a 
skyscraper, or protrusion, effect, which is even more pronounced for lower diffusion rates, as 
discussed in Appendix |A] Since the protrusions contain a significant amount of molecular 
hydrogen and since we have not seen them up to now, it is likely that they can only be asso- 
ciated with cold dense cloud conditions with the diffusion to desorption of 0.5 as is used here. 
The width and height of the protrusions are associated with the random angle deposition 
method. With this method, the protrusions continue to grow since they prevent incoming 
atoms from reaching the surface via a shadowing effect. It is likely that the tendency to form 
such noticeable skyscraper structures will be mitigated in the interstellar medium for several 
reasons. First, the grain surfaces are curved, which reduces the shadowing effect. A standard 
grain with a radius of 0.1 /im has a circumference of approximately 2000 sites, which means 
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the snapshots in Figure [5] represent arcs of 9°. For a small grain of r = 0.02/mi, it corre- 
sponds to 45°. Secondly, the rate of diffusion considered in our models is likely a lower limit 
because diffusion also occurs in the bulk structures that are formed. These will result in a 
smoother surface with less molecular hydrogen (see the following paragraph). To confirm our 
hypothesis concerning the random angle deposition method, we also performed simulations 
with perpendicular deposition. Large single protrusions containing a high percentage of H 2 
were not obtained under these conditions; only structures similar to the left panel in Fig. [5] 
were found. The sticking fraction (fraction remaining in the mantles) in the simulation with 
the perpendicular deposition is much higher and a thicker layer of ice is formed. Angular 
deposition gives more irregular structures at the surface that have a weaker binding energy 
and can therefore desorb more easily, resulting in a lower sticking fraction. 



Our results can profitably be compared with another simulation. iKimmel et al.l (120011 ) 
studied the porosity of ice mantles formed by deposition of water molecules using a bal- 
listic algorithm. In their algorithm they deposited the molecules under specific angles to 
simulate a molecular beam or with varying angles to simulate a background pressure. One 
of their conclusions is that ices deposited from a background pressure are very rough and 
porous, in agreement with our result that random angle deposition results in such structures. 
Kimmel et al.l (120011 ) also let a molecule anneal immediately after its deposition. Their an- 
nealing algorithm consists of several cycles in which the deposited species hop to energetically 
more favored sites. Annealing results in denser ices (fewer narrow skyscrapers). We do not 
include annealing in our simulations, but we do include diffusion of H and O atoms whereas 
the water molecules of IKimmel et al.l (120011) do not diffuse thermally. In agreement with 



the observations of IKimmel et al.l (120011 ) is our result that more diffusion results in denser 
ices, as can be seen by comparing the broad skyscrapers in Figure [5] with the narrow ones 
discussed in Appendix |A] 



5. DISCUSSION 

Using the CTRW Monte Carlo method, we have shown that mantles of water and 
other ices can grow efficiently in cold and dense regions of the interstellar medium whereas 
only a small amount, in the sub monolayer regime, is formed under the conditions of diffuse 
interstellar clouds. Ices in translucent sources occupy an intermediate regime. We have also 
looked carefully at the morphology of the ices. This initial exploration is, however, dependent 
on a variety of poorly constrained parameters such as barriers to diffusion, which are generally 
less well understood than adsorption energies. With our standard model, the ices appear 
to be much denser in diffuse and translucent sources, although there are many pores in the 
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interior. Also, there is a vertical stratification in these ices, in which water-ice dominates in 
the uppermost monolayers while unusual species such as OH and O3 are abundant in the lower 
monolayers. In the dense sources, the ices, which are dominated by water and occasionally 
H 2 , appear to develop a skyscraper- type structure, which is somewhat dependent on the rate 
of diffusion adopted. This very rough structure is probably exaggerated somewhat in our 
calculations because we do not properly include diffusion inside the bulk. Although the water 
ice in the diffuse and translucent clouds is produced by the association reaction between H 
atoms and OH radicals, which is customarily assumed to be the dominant source of water, the 
situation in the dense clouds, which contain low abundances of gas-phase hydrogen atoms, 
is quite different. 

Observations in dense cloud assemblies such as Taurus show no water ice below a thresh- 
ol d value of Ay = 3. 2 and linear growth of the column density of water ice above this value 



(IWhittet et al.ll200ll ). In particular, with a least-squares- fit, the expression 

iV(H 2 o) = g« s - A th ) (15) 
(1.30 ± 0.04 x 10 17 ) cm" 2 and A th = 



can be determined (IWhittet et al.ll200ll ). where q 
3.2 ± 0.1, using the relationship 



iV(H 2 0) 



U 



x 10 18 [cm- 2 ]r 3 .o. 



(16) 



Relating these observations to the results of our simulations is not completely straightfor- 
ward. The measured threshold value of Ay = 3.2 is observed through the cloud and gives an 
edge-to-edge value. The observed ice content is the integrated contribution along a line-of- 
sight that typically contains heterogeneous material. The visual extinction experienced by 
most material along this line is lower, on average one-half of the edge-to-edge value. 

We make the assumption that the extinction along any line of sight is dominated by 
one object, defined by the physical conditions used for our models and listed in Table HI 
The number of monolayers at 10 5 yr for the dense objects and 3 x 10 5 yr for the diffuse 
and translucent objects is converted into column den sities of ice us i ng th e gas-to-dust ratio 
for standard-sized grains and typical cloud depths (jNguyen et all 120021 ) . The results are 
plotted in Figure [6] at the six different visual extinctions used in our simulations. The 
several different results at each extinction correspond to the different physical conditions 
chosen. Also depicted in the figure are two dashed lines that represent the results of equation 
(IT5|) assuming that the observed edge-to-edge extinctions correspond to our extinctions - 
Ay m = Ay s - and that they correspond to twice our extinctions: Ay m = \Ay s . Finally, 
there is a solid line representing roughly the lower detectable limit for water ice based on the 
H 2 optical depth at 3/m towards Cygn OB 2 no. 12 of r < 0.02 Jwhittet et alJll997h . This 



optical depth corresponds to an iV(H 2 o) of 3.6 x 10 16 cm 2 , via Eq. [TBI which is on the order 
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of a few monolayers. Assuming that the typical extinction seen by ice along a line-of-sight is 
indeed 1/2 the edge-to-edge value, our results should approach the upper dashed line shown 
in Figure EJ The plot indicates, however, that our simulations produce less ice than needed 
by a factor of ~ 3 — 10. 

This shortfall can be accounted for by increasing the depth of our sources and, for 
dense clouds, extending the simulations to later times since the calculated ice abundances 
are still increasing rapidly at our latest time of 10 5 yr. However, cosmic ray desorption 
start s to become important after 10 5 yr and part of the grain mantle can desorb in this 
way ( Herbst &: Cupper] 2006 ). It may be that we underestimate the efficiency of water-ice 



formation, especially at intermediate values of the visual extinction. If so, we can produce 
more ice if the rates of diffusion are enhanced, as discussed in Appendix |A] Also a decrease 
of the shadowing effect due to the presence of curvature, as discussed earlier, can result 
in thicker ice layers on small grains as we saw in the perpendicular deposition simulations. 
Yet another possibility is that some ice is indeed accreted from gaseous water, possibly in 
the cold gas following a shock, which would allow water to be produced in the gas via the 
reaction between OH and H2, despite its activation energy. 

E.H. would like to thank the National Science Foundation (US) for supporting his re- 
search program in astrochemistry. 
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A. THE INFLUENCE OF THE HOPPING BARRIERS 



A ratio of 0.5 between the diffusion barrier and desorption energy has been assumed 
throughout the paper. This appendix sho ws the simulation r esults using an alternative ratio 
of 0.78, generalizing the result obtained by lKatz et al.l (119991 ) for atomic hydrogen on olivine. 



The results are shown in Figures [7] to [10] and in Table [71 Figure [7] shows surface 
abundances for diffuse and translucent cloud conditions as a function of time, while Figure 
[8] does the same for dense cloud conditions. Figure M contains vertical cross sections of ice 
under diffuse and translucent conditions, while Figure [TD] contains such images for dense cloud 



This preprint was prepared with the AAS IATgX macros v5.2. 



-21 - 



conditions. Table [7J lists H2 recombination efficiencies for diffuse and translucent clouds. As 
can be seen by comparing these figures and table with the analogous ones for our standard 
hopping barrier, there are important differences in the amount of ice, the recombination 
efficiency for molecular hydrogen formation, and in the structure of the ice between these 
results and the previous ones. Let us consider them in turn. 

Considering the structure of the mantles, the ice is less dense for the higher hopping 
barrier, as can be clearly seen by comparing Figures [3] and M and Figures and [TDl This 
can be explained by realizing that for the lower hopping barrier, oxygen atoms can diffuse 
to more favorable positions (with more neighbors) before reaction. 

Large protrusions as the one shown in the left panel of Figure [10] can be observed for 
both diffusion ratios, but they are more common and higher for the higher ratio, and are 
even present in translucent sources. The main reason is that the mobility of species like 
O, H, and H2 is limited to such an extent that they cannot move to energetically more 
favorable position before it is covered by another particle. In cold dense clouds this effect 
is the strongest since the surface temperature is low and the flux of, especially, molecular 
hydrogen is very high. As explained in the main text, the formation of these protrusions is 
amplified by self-shadowing of the surface by the protrusion. Since a real interstellar grain 
surface will have a curvature, which becomes significant for smaller grains of approximately 
0.02 fim, this shadowing effect will be less in reality. 

More ice (see Figure [7J in comparison with Figure [2]) and more molecular hydrogen are 
formed for the higher ratio. This is consistent with previous results for molecular hydrogen 
formation. If we compare the right columns in Tables [5] and [7J we see that also in the case 
where only a hydrogen flux is present, the production of molecular hydrogen is much higher. 

B. THE INFLUENCE OF THE ROUGHNESS OF THE GRAIN 

This appendix contains the results of simulations starting with a less rough grain than 
our standard one. With this smoother surface, the running time of the simulations becomes 
significantly longer. This effect occurs because on the smoother surface the average hopping 
time is smaller so that the simulation program spends more of its time evaluating hopping 
events. For this reason we performed simulations only for a few of the conditions in Table 
SJ Surface abundances vs time are shown in Figure [TT] for panels A3, B3, C3, and D3. 

In general, there is little difference between these results and our standard ones. The 
amount of ice that is formed is very similar for all simulations, except for the initial peak 
around 5 x 10 4 yr for B3 when the smoother surface is used. The surface abundances of the 
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other species change slightly in the very diffuse regions. A3 has less H 2 and B3 less O for the 
smoother surface. In the other two cases, the surface is completely covered with ice, which 
determines the energetics so that the difference in surface abundances between the smoother 
and the rougher surfaces becomes negligible. 




Fig. 1. — The surface topology of the bare carbonaceous grain, r epresented by a lattice o f 
adsorption sites. The lighter colors indicate the higher areas. See ICuppen fc Herbstl (120051 ) 
for a detailed discussion of how the surface was generated. 




Fig. 2. — The surface abundance as a function of time for different conditions in a diffuse or 
translucent molecular cloud. See Table H] for physical conditions in each panel. 
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Fig. 3. — Vertical cross sections of the ice mantles after 3 x 10 5 yr for four different physical 
conditions in diffuse and translucent clouds (A2, B2, C2, D2; Table HJ). The molecules are 
color coded. 
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Fig. 4. — The surface abundance as a function of time for different 
molecular cloud. The notation a(b) implies a x 10 6 . See Table H] for 
each panel. 



conditions in a dense 
physical conditions in 
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Fig. 5. — Vertical cross sections of the ice mantles after 1 x 10 5 yr for two different dense-cloud 
physical conditions. Left panel (E2): A v = 5, T gr = 12 K, T gas = 20 K and n H = 5 x 10 3 
cm" 3 . Right panel (F2): A v = 10, T gr = 10 K, T gas = 10 K and n u = 2 x 10 4 cm" 3 . The 
molecules are color-coded. 
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Fig. 6. — Column density of water ice at the end of the simulations (3 x 10 5 yr for A\ 
0.5, . . . , 3 and 10 5 yr for Ay m = 5 and 10) versus the extinction used in our simulations for 
all physical conditions given in Figs. 
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Fig. 7. — The surface abundance with a high hopping barrier-to-desorption energy of 
as a function of time for different conditions in a diffuse or translucent molecular cloud. 



-30- 




Fig. 8. — The surface abundance with a high hopping barrier-to-desorption energy of 0.78 
as a function of time for different conditions in a dense molecular cloud. See Table H] for 
physical conditions in each panel. 
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Fig. 9. — Vertical cross sections of ice mantles after 3 x 10 5 yr for four different physical 
conditions in diffuse and translucent clouds (A2, B2, C2, D2; Table H]) obtained with a high 
hopping barrier-to-desorption energy of 0.78. The molecules are color-coded. 
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Fig. 10. — Vertical cross sections of ice mantles after 1 x 10 5 yr for two different dense-cloud 
physical conditions obtained with a high hopping barrier-to-desorption energy of 0.78. Left 
panel (E2): A v = 5, T gr = 12 K, T gas = 20 K and n H = 5 x 10 3 cm" 3 . Right panel (F2): 
Ay = 10, T gr = 10 K, T gas = 10 K and «h = 2x 10 4 cm -3 . The molecules are color coded. 
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time (yr) 

Fig. 11. — The surface abundance as a function of time starting with a bare grain with a 
smooth surface structure. 



-34- 



Table 1. SURFACE REACTIONS 



Surface reactions 





Reaction 




& 


E a (K) 


H - 


f H 


-» H 2 




0.991 a 





H - 


h O 


-» OH 




0.991 a 





H - 


f- OH 


-* H 2 




0.991 a 





O - 


t- O 


-* O2 




0.991 a 





H - 


t- 2 


•4 2 H 




0.991 a 


1200 b 


H - 


f 2 H - 


-» H2O2 




0.991 a 





H - 


f H2O2 - 


-+ H 2 


+• OH 


1 


1400 c 


H - 


t- O3 


-» 2 


4- OH 


1 


450 d 


H 2 - 


h OH 


-» H2O 


+ H 


1 


2600 c 


O - 


f 2 


- O3 




0.991 a 






a Based on photod i ssocia tion of water results by 
iKroes fc Anderssonl j2006h . 

iMelius fc Blind Jl979h 

telemm et all Jl975h 

^Lee et alj Jl978^ 

jSchifj dl973h 
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Table 2. PHOTODISSOCIATION REACTIONS 



Reaction 51 






&p hoto (s 


) 'Yphoto 


a cr 


OH -> 


O - 


h H 


1.68(-10) 


1.66 


1.02(3) 


H 2 -> 


H - 


h OH 


3.28(-10) 


1.63 


1.94(3) 


O2 -> 


O - 


h O 


3.30(-10) 


1.4 


1.50(3) 


2 H -> 


O - 


h OH 






1.50(3) 


O2H -> 


H - 


h O2 






1.50(3) 


H 2 2 -> 


OH - 


h OH 






3.00(3) 


O3 -» 


O2 - 


1- 


1.9(-9) 


1.85 b 





a gas-phase photodissociation rates and products are 
used 

Ivan Dishoeck et al.l [|2006ri 
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Table 3. DESORPTION ENERGIES IN KELVINS 



Substrate 



Absorbate 


carbon 


H 


O 


OH 


H 2 


O2 


H2O 


O3 


O2H 


H2O2 


H 


658 








45 


45 


650 


45 




45 


O 


800 






55 


55 


55 


800 


55 


55 


55 


OH 


1360 




240 


240 


240 


240 


3500 


240 


240 


240 


H 2 


542 


30 


30 


30 


23 


30 


440 


30 


30 


30 


Oa 


1440 


69 


69 


69 


69 


69 


1000 


69 


69 


109 


H 2 


2000 


390 


390 


390 


390 


390 


5640 


390 


390 


390 


O3 


2240 


120 


120 


120 


120 


120 


1800 


120 


120 


120 


O2H 


2160 




300 


300 


300 


300 


4300 


300 


300 


300 


H 2 2 


2818 


340 


340 


340 


340 


340 


4950 


340 


340 


340 



Note. — See discussion in text 
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Table 4. PHYSICAL CONDITIONS OF CHOSEN MODELS 



Panel 


Ay 




T gas (K) 


nn (cm 3 ) 


Form of H 


Al 


0.5 


20 


100 


1.0(2) 


11 


A2 


0.5 


18 


80 


1.0(2) 


H 


A3 


0.5 


16 


60 


1.0(2) 


H 


A4 


0.5 


18 


80 


5.0(2) 


H 


Bl 


1.0 


18 


80 


2.5(2) 


H 


B2 


1.0 


16 


60 


2.5(2) 


H 


B3 


1.0 


14 


40 


2.5(2) 


H 


B4 


1.0 


16 


60 


5.0(2) 


H 


CI 


2.0 


17 


70 


5.0(2) 


H 


C2 


2.0 


15 


50 


5.0(2) 


H 


C3 


2.0 


13 


30 


5.0(2) 


H 


C4 


2.0 


15 


50 


1.0(3) 


H 


Dl 


3.0 


16 


60 


1.0(3) 


H 


D2 


3.0 


14 


40 


1.0(3) 


H 


D3 


3.0 


12 


20 


1.0(3) 


H 


D4 


3.0 


14 


40 


1.5(3) 


H 


El 


5.0 


14 


40 


5.0(3) 


H 2 


E2 


5.0 


12 


20 


5.0(3) 


H 2 


E3 


5.0 


10 


10 


5.0(3) 


H 2 


E4 


5.0 


12 


20 


1.0(4) 


H 2 


Fl 


10.0 


12 


20 


2.0(3) 


H 2 


F2 


10.0 


10 


10 


2.0(4) 


H 2 


F3 


10.0 


10 


10 


1.0(4) 


H 2 


F4 


10.0 


10 


10 


5.0(4) 


H 2 



Note. - 



The notation a(b) implies a X 10 6 . 
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Table 5. H 2 RECOMBINATION EFFICIENCY 



12 3 4 



A 


0.05 


0.02 


0.04 


0.17 


0.22 


0.57 


0.05 


0.11 


B 


0.06 


0.28 


0.15 


0.59 


0.24 


0.72 


0.17 


0.60 


C 


0.24 


0.54 


0.29 


0.65 


0.48 


0.71 


0.22 


0.68 


D 


0.11 


0.61 


0.43 


0.74 


0.52 


0.71 


0.54 


0.73 



Note. — The left columns arc results in the presence 
of an oxygen flux, the right in its absence. 
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Table 6. CONTRIBUTIONS OF THREE REACTIONS TO WATER FORMATION 







1 






2 






3 






4 




A 


99.9 


0.0 


0.1 


99.6 


0.0 


0.4 


95.9 


0.0 


4.2 


99.5 


0.0 


0.5 


B 


99.5 


0.0 


0.5 


98.3 


0.0 


1.7 


96.4 


0.0 


3.6 


98.5 


0.1 


1.4 


C 


98.4 


0.5 


1.1 


97.1 


0.8 


2.0 


93.1 


1.0 


5.9 


95.9 


1.3 


2.9 


D 


97.6 


1.4 


1.0 


93.4 


1.4 


5.2 


90.5 


1.6 


7.9 


93.2 


0.9 


5.9 


E 


14.3 


21.0 


64.7 


11.6 


19.3 


69.1 


7.0 


16.7 


76.3 


9.7 


18.3 


72.0 


F 


10.1 


18.7 


71.2 


6.5 


16.7 


76.8 


6.2 


16.5 


77.2 


6.7 


16.2 


77.1 



Note. — From left to right: H + OH -<• H 2 0, H + H 2 2 ->■ H 2 + OH, 
and H 2 + OH ->■ H 2 + H. 
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Table 7. H 2 RECOMBINATION EFFICIENCY WITH SLOW DIFFUSION 



12 3 4 



A 


0.47 


0.56 


0.58 


0.72 


0.61 


0.68 


0.59 


0.69 


B 


0.67 


0.74 


0.64 


0.67 


0.68 


0.71 


0.67 


0.68 


C 


0.40 


0.71 


0.43 


0.74 


0.49 


0.65 


0.37 


0.75 


D 


0.30 


0.70 


0.34 


0.68 


0.32 


0.69 


0.31 


0.67 



Note. — The left columns arc results in the presence 
of an oxygen flux, the right in its absence. The calcu- 
lations are done with a ratio of 0.78 between hopping 
barrier and desorption energy. 



